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ABSTRACT 

We present a 1.1mm emission map of the OMC1 region observed with AzTEC, a new large-format 
array composed of 144 silicon-nitride micromesh bolometers, that was in use at the James Clerk Maxwell 
Telescope (JCMT). These AzTEC observations reveal dozens of cloud cores and a tail of filaments in a 
manner that is almost identical to the submillimeter continuum emission of the entire OMC1 region at 
450 and 850 fim. We perform Fourier analysis of the image with a modified periodogram and the density 
power spectrum, which provides the distribution of the length scale of the structures, is determined. 
The expected value of the periodogram converges to the resulting power spectrum in the mean squared 
sense. The present analysis reveals that the power spectrum steepens at relatively smaller scales. At 
larger scales, the spectrum flattens and the power law becomes shallower. The power spectra of the 
1.1mm emission show clear deviations from a single power law. We find that at least three components 
of power law might be fitted to the calculated power spectrum of the 1.1mm emission. The slope of 
the best fit power law, 7 « —2.7 is similar to those values found in numerical simulations. The effects 
of beam size and the noise spectrum on the shape and slope of the power spectrum are also included 
in the present analysis. The slope of the power law changes significantly at higher spatial frequency as 
the beam size increases. 
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1. INTRODUCTION 

The Orion Nebula is located in the northern part 
of the Orion A molecular cloud and consists of OMC1, 
OMC2, OMC3, and OMC1-S. These components are 
smoothly connected and form integral shaped filaments 
seen in continuum observations at sub-millimeter, in- 
frared, and optical wavelengths. The M42, Orion Neb- 
ula, is known to be the best studied star-forming region 
in the nearest giant molecular cloud (GMC) where mas- 
sive stars have formed and has been extensively studied 
at all wavelengths. OMC1 in the Orion A molecular 
cloud possesses the most massive component contain- 
ing at least 2,000 young stellar objects (YSOs) (O'Dell 
2001) and is strongly connected to the surrounding gas. 
There is also a report on the VIA NH3 observations 
and revelations of the filamentary and clumpy struc- 
ture of OMC1 in the Orion A molecular cloud (Wise- 
man & Ho 1996, 1998). OMC1 has an on-going mas- 
sive star formation and strong interactions between gas 
and young stellar objects (YSOs), their bipolar stellar 
jets, and outflows (Mann & Williams 2010; Bally et 
al. 2011). As we have mentioned above, the complex 
and filamentary structure of the Orion A molecular 
cloud, especially the Orion Integral Shaped Filament 
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was revealed by sub-millimeter continuum observations 
by Lis et al. (1998) and Johnstone & Bally (1999). 
The condensation mass spectrum was also explored and 
temperature distribution was examined (Johnstone & 
Bally 1999). Mookerjea et al. (2000) found that the 
coldest clump in the Orion A molecular cloud had a 
temperature of about 15 K. These cold clumps can 
be traced at the continuum observations using the ob- 
servations at mm wavelength. Despite various obser- 
vations of the Orion A molecular cloud over a wide 
range of wavelengths for over a decade including the 
SEST SIMBA (Nyman et al. 2001), thermal contin- 
uum emission at 1.1mm has been observed for the first 
time only recently, owing to the development of the As- 
tronomical Thermal Emission Camera (AzTEC). This 
AzTEC camera is composed of an array of 144 nitride 
micro-mesh composite bolometers (Wilson et al. 2008). 
This camera was originally designed for a millimeter- 
wavelength bolometer camera to be installed on the 
Large Millimeter Telescope (LMT). But in year 2005 
and 2006, it was installed on the James Clerk Maxwell 
Telescope (JCMT) in Hawaii, Mauna Kea and per- 
formed biased and unbiased surveys, of the northern 
sky. 

The present study analyzes the 1.1mm AzTEC ob- 
servations of the OMC1 components and reports the 
fact that the 1.1mm emission from OMC1 is distributed 
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in a manner almost similar to that of the 850 /im emis- 
sion (Johnstone & Bally 1999) observed with the Sub- 
millimeter Common-User Bolometer Array (SCUBA) 
on the JCMT. In the present analysis, we performed 
a power spectrum analysis of the 1.1mm intensity map 
which is dominated by thermal emission from cold dust 
in OMC1 in the Orion A molecular filaments. Since the 
dust grains are coupled with gas, elucidating the struc- 
ture of the dust emission will also play a crucial role in 
understanding the dynamics and the structure of the 
gas. In general, the density fluctuations of the inter- 
stellar dust (Draine & Lazarian 1998) occur in the cold 
and dense interstellar medium (ISM). 

The power spectrum analysis using the Fourier con- 
version of the given images is especially important be- 
cause they can reveal the statistical properties of struc- 
tures which may be present across the images. The gen- 
eration of the hierarchical structures or fractal struc- 
tures (Falgarone et al. 1991; Elmegreen & Falgarone 
1996) often presents across a range of spatial scales in 
the ISM. It has also been suggested that this could 
be attributed to turbulent emission since this can in- 
crease gas density and initiate star forming instabilities 
(Mac Low 2004; Burkert 2006). However, the large 
scale properties of turbulence in the ISM, including 
those of the gas and dust grains are still poorly un- 
derstood. But there has been an explosion of research 
in this area. In general, the energy input from star for- 
mation is thought to be a major driving force for the 
turbulence in the ISM (Scalo & Elmegreen 2004). It 
is also well known that measurements of the slopes of 
the power spectra (Crovisier & Dickey 1983; Gautier 
ct al. 1992) from the intensity fluctuations in the tur- 
bulent medium can provide an insight into the nature 
of scales present in the hierarchical or non-hierarchical 
structures (Blitz & Williams 1997; Hartmann 2002) in 
the interstellar medium (Brunt 2010). Thus, the re- 
sults from the present studies of the power spectrum 
analysis can reveal the nature of the structures seen in 
the 1.1mm image of OMC1 in the Orion A molecular 
cloud and filaments. The results of this present analy- 
sis are presented in Section 3 and 4. Section 2 provides 
procedures of image processing and the resultant data. 

2. OBSERVATION AND DATA PROCESS- 
ING 

In this section, we describe the processing of the 
1.1mm AzTEC survey data sets. The observations of 
OMC1 were performed in raster-scan mode between 28 
December 2005 and 1 January 2006 using the AzTEC 
mounted on the JCMT, giving a beam size of 18 arcsec 
and a scan velocity of 120 arcsec s _1 . A total inte- 
gration time of 102.9 seconds was used. From a to- 
tal of ten scan maps, five scans were combined, cover- 
ing an observing area of 655.35 arcmin 2 , which corre- 
sponds to about 3.14 x 3.14 pc 2 . Average RMS was 
measured to be 0.22 mJy. The observed region was 
centered on RA=05 h 35 m 14.4 s , DEC=-05°22'32.3" in 



J2000. AzTEC suffered a pointing error of 2". Beam 
maps were made from observations of Uranus with a 
mean flux density of 42.2 Jy at 1.1mm during the 
JCMT observing run (Wilson et al. 2008). 

A pointing model was generated using CRL618, a 
post-asymptotic giant branch (post-AGB) star with a 
mean flux density of 2.7 Jy at 1.1mm. The inten- 
sity calibration errors were estimated to be about ±6- 
13% (Wilson et al. 2008). The raw bolometer data 
was calibrated and edited in the IDL AzTEC software 
(Perera et al. 2008; Scott et al. 2008; Wilson et al. 
2008; Austermann et al. 2009). Each data set was re- 
garded as being in the time sector using the data stor- 
age method. Since all the data was stored as a func- 
tion of time, they were also referred to as time-stream 
data. The AzTEC data set was also in the spatial sec- 
tor. Each observed data set i, Mi, was composed of 
the atmospheric signal, A,, the astronomical signal, S, 
and the noise, Ni, in the spatial sector. An astronom- 
ical signal estimator, Si, was built by cleaning the raw 
data rrii, and the time-domain analogs of the Mi. 

For extended sources, 5* correlates with the time- 
stream signals as Aj. The atmosphere a>i, time-domain 
analogs of Aj in the time sector, may not be orthogonal 
to the true sky signal in the time sector, s. An iterative 
technique was developed to maximize the orthogonality 
between the principal components being cut from time- 
stream data and the real astronomical signal. This 
means that measurement of fluxes in the initial map 
was subtracted from the raw data set. Then, the data 
streams were cleaned and mapped again. Any signif- 
icant detections in the second map were added to the 
source model of the first generated map. The summed 
model was then subtracted from the initial bolometer 
signals. In summary, an iterative cleaning algorithm 
was performed using the following procedure: 1) set of 
Si was found from the set of ra,. S was created by 
co-adding all the Si; 2) f,, the estimator for the time- 
stream data unrelated to the astronomical signal, was 
created by subtracting s from m,; 3) ?{ was cleaned 
and maps were created for the set of Rf, these were 
co-added to find R, the residual map of the sky signal; 
4) a new estimator of the true sky signal, S' was cre- 
ated by adding the residual map, R to the map from 
step 1); 5) we returned to the step 2) and subtracted s' 
from the input time streams. The final AzTEC image 
of OMC1 was the result of thirty iterations of clean- 
ing as described above. Any induced artifacts from the 
iterative cleaning processes were compensated for by 
the inverse effects of any subtracted non-real features 
in the time stream data. This process was repeated 
thirty times until the result converged to the model- 
subtracted map. 
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5 1, 3o m 30 ! 35™20 s 35 m 10 s 

Right Ascension (J2000) 

Fig. 1. — The AzTEC 1.1mm emission from OMC1 in the 
Orion A molecular cloud. The peak flux of the 1.1mm emis- 
sion at Orion KL of OMC1 in the Orion A molecular cloud 
is about 63.9 Jy/beam. The calibration error is reported as 
6 - 13% of the flux in the paper by Wilson et al. (2008). 



3. FOURIER ANALYSIS 

We present observations of the 1.1mm emission from 
OMC1 in Fig. 1. The thermal dust emission from cool 
dust at 1.1mm can be the result of heating from in- 
termediate mass stars (5—20 M©) and the interstellar 
radiation field. The morphology of the emission is sim- 
ilar to that of the 850 /im emission (Johnstone & Bally 
1999) seen with SCUBA at JCMT. The northern peak 
of OMC1 is associated with Orion KL (Kleinmann & 
Low 1967; Allen et al. 1993; Menten et al. 1995; Lee 
& Burton 2000; Zapata et al. 2011) in OMC1 and the 
southern peak of OMC1 corresponds to the Orion S 
source (Schmid-Burgk et al. 1990). The peak fluxes of 
these two sources in Fig. 1 were measured in Jy/beam 
and presented in Table 1. 

To extract quantitative information from the sta- 
tistical properties of the interstellar emission, we per- 
formed power spectrum analysis of the intensity map 
shown in Fig. 1. This analysis provides a useful infor- 
mation on understanding the physics and structure of 
the interstellar dust medium which was dominated by 
thermal emission from the cold dust in OMC1 and fil- 
aments. The parameters of the images from the power 
spectrum analysis were often used to study the sta- 
tistical properties of pre-stellar cores in the molecu- 
lar clouds and filaments as well as in the cosmic mi- 
crowave background (CMB) (Shaw & Lewis 2012) and 
cosmic infrared background (CIB) radiation analysis 
(Finkbeiner et al. 1999; Keisler et al. 2011). The 
power spectrum has also been used extensively in struc- 



tural analysis of the images (Combes et al. 2012; see 
references therein). The main goal of this analysis is 
to place some constraints on the density structure of 
the interstellar dust emission arising from OMC1 using 
Fourier analysis (Bracewell 1986). Since dust grains 
are coupled with gas, elucidating the structure of the 
dust emission will also play an important role in un- 
derstanding the dynamics and the structure of the gas. 
The power spectrum of the 1.1mm dust emission from 
OMC1 is shown in Fig. 2. 

The power spectrum, -P(fc), was measured by carry- 
ing out the Fourier transform (Bracewell 1986; Press 
et al. 1992) of the auto-correlation of the image (Kim 
et al. 1999). The Fourier transform of the product 
of the data with a window function reflects the con- 
volution of the Fourier transform of the data with the 
window function (Bracewell 1986; Press et al. 1992). 
In the present analysis, a nonparametric periodogram 
estimator was used in the Fourier transform (Eq. 1) of 
the image (Fig. 1). 




Fig. 2. — Here, the density power spectrum of the AzTEC 
1.1mm emission from OMC1 is presented. The estimate of 
the power spectrum was calculated with different aperture 
diameters of 18" (blue), 36" (green solid), 56" (red solid), 
and the noise (black) for the comparison. We also compared 
the results of the analysis for the beam convolution with 
those for the averaging of the pixels (dashed lines) in the 
image. 

The periodogram estimate of the power spectrum, 
P(w, N) is defined on N/2 + 1 frequencies in Eq. 1 



P(w,N) = 



N-l 

£ 

k=-N+ 



(i) 
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Table 1. 

Emission properties of Orion KL and Orion S of OMC1 in the Orion A Molecular Cloud at 1.1mm 





AzTEC Coordinate 


Peak Flux 




RA[J2000] DEC[J2000] 


[Jy beam -1 ] 


Orion KL 


05:35:14.17-05:22:33.35 


63.9 


Orion S 


05:36:13.52 -05:24:07.35 


26.8 



where 

N-l-k 

r *^=N £ <n + k)x*{n). (2) 

71=0 

The data were partitioned into K segments and each 
2M consecutive points at a number of discrete frequen- 
cies. Here, N is 2M. k refers to a wavenumber. The pe- 
riodogram estimates were averaged at each frequency, 
which reduced with the variance of the estimates by 
a factor of K. Welch's modified periodogram is justi- 
fied in Eq. 3. The expected value of the periodogram 
converges to the power spectrum in the mean squared 
sense. 



PH « ^ E \ J2 w ( n >( n + iD ) e ~ 2 ™ w \ 2 ( 3 ) 

i=0 71=0 

where the overlap occurs at L — D data points and 
N = L + D(K - 1) (Press et al. 1992). w(n) refers to 
a window. 

The periodogram was calculated by computing the 
discrete Fourier transformation. Here, the result was 
an array of power measurements for each frequency bin. 
The amplitude of each frequency bin was determined by 
centering a window for an each bin and measuring the 
amount of signal falls within the window. The individ- 
ual data sets are commonly overlapped in Welch's peri- 
odogram and this windowing of the segments makes a 
modified periodogram. This modified periodogram in- 
cluding overlap processing computes the weighted aver- 
age of the sum of the spectrum from each segment with 
the window. The window functions usually have more 
influence on the data at the center of the set than on 
the data at the edges. In calculating the periodogram, 
50% of overlapping is treated in this paper. 

The variance of the periodogram estimate can be 
defined as 

a 2 = i |>(«)- Ml 2 (4) 

71=0 

where \i refers to a mean. When the segments are 
overlapped by one half of their lengths, the variance 



of Welch's periodogram reaches an optimal value. This 
reduces the variance by a factor of about 0.818-KT (Press 
et al. 1992). If there is an increase in the number of 
subsequences to be averaged, then the variance of the 
estimate is reduced significantly. Therefore, we can say 
that the power spectrum estimation is approximately 
unbiased. 

4. POWER SPECTRA 

The power spectrum of density fluctuations can be 
related to the power spectrum of intensity fluctuations 
(Bondi et al. 1994). In this analysis, we found that it 
was difficult to fit a single power law to the observed 
spectrum in Fig. 2. The power law at large scales is 
shallow because it connects to the large-scale portion, 
while the power spectrum is dominated by the size and 
shape of the map. The power spectrum obtained from 
the intensity map of OMC1 in the present analysis dis- 
plays at least three distinct power law exponents as 
can be seen in Fig. 2. The best-fit slope in the range 
30 £ k £ 100 is 7 w -2.66 ± 0.3. This result is com- 
parable to the spectral index 7 w —2.7 that was found 
in numerical simulations (Padoan et al. 2004) within 
error estimates. There are distinct spectral breaks at 
spatial frequencies of k s=s 30 and at k w 100 possibly. 
While the break in the power spectrum occurring at the 
thickness of the galactic disk was noted in Elmegreen, 
Kim, & Staveley-Smith (2001), the preferred scale ap- 
parent in the power spectrum was also discussed by 
Brunt (2010) and in the previous studies by Blitz and 
Williams (1997) and Hartmann (2002). 

At smaller scales the observed power spectrum of the 
image representing the 1.1mm emission from OMC1 
approaches the noise level. The power spectrum steep- 
ens at relatively smaller scales. The power spectrum 
approaches the noise level at higher spatial frequency 
(Fig. 2). Estimates of the power spectra at smaller 
scales can often underestimate the value, since struc- 
tures can be suppressed below the sonic scale (Vazqucz- 
Semadeni, Ballcstcros-Paredes, & Klessen 2003). From 
the present analysis (Fig. 2), the power spectrum 
steepens at relatively smaller scales. The logarithmic 
slopes of the power-law spectra obtained at low spa- 
tial frequency seem to be shallower than the logarith- 
mic slope of Kolmogorov spectrum. A Kolmogorov 
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spectrum has a power law with a logarithmic slope of 
— 11/3, conventionally taken as evidence for turbulent 
fluctuations in dynamical energy exchange among ir- 
regularities of different scales. The spectral indices of 
the power spectra of the dust emission derived from 
other studies (Gautier et al. 1992) are also rather shal- 
lower than that of Kolmogorov spectrum. Most spec- 
tral indices measured in other studies such as Cirrus 
cloud studies and Cassiopeia A (Roy et al. 2008) range 
from —3.6 to —2.2. But in the present Fourier anal- 
ysis, the slope of the power law at high spatial fre- 
quency becomes steeper than the logarithmic slope of 
Kolmogorov spectrum. Often spectra with logarithmic 
slope of four or more are considered to be caused by 
deterministic structures. However, it is difficult to ex- 
plain the power source for forming these substructures. 

Recent studies of turbulence indicate a rather shal- 
low slope of the density power spectrum. For example, 
Kritsuk & Norman (2004), Beresnyak, Lazarian, & Cho 
(2005), Kim & Ryu (2005), and Bournaud et al. (2010) 
found slopes which were rather smaller than 2 from 
their magnetohydrodynamic (MHD) turbulence. The 
slopes of power law become smaller as the Mach num- 
ber increases in their simulations and measurements. 
In the present Fourier analysis, the slope of the power 
law at low spatial frequency, k <; 30 is 7 w —1.5 ± 0.1 
and smaller than 2. In the transonic turbulence model, 
the density distribution is often generated by disconti- 
nuities created by weak shocks overlaid on the turbu- 
lent background. In the supersonic turbulence, the den- 
sity distribution is often characterized as peaks, or mass 
concentrations generated by strong shocks (Williams, 
Blitz, & McKee 2000). As the Mach number increases, 
the power spectrum flattens (Kim & Ryu 2005), and 
density concentrations appear as sheets and filaments 
(McKee & Ostriker 2007). To interpret these results in 
terms of the density structure (Brunt & Heyer 2002) 
and in comparison with deductions from gas tracers, 
we reproduced the results by considering a model spec- 
trum based on the 1.1mm emission which is dominated 
by the gray-body emission from large dust grains at 
thermal equilibrium temperature (Desert et al. 1990). 
Here, the large grain equilibrium temperature is re- 
lated to the local radiation field strength (Stepnik et 
al. 2003) and spectrum which depends on the presence 
or absence of nearby heat sources as well as extinc- 
tion. Since grain structure variation and size can affect 
the emissivity coefficient and can result in a variation 
of the equilibrium dust temperature, this will conse- 
quently constrain the pattern of the power spectra at 
smaller and intermediate scales in the interstellar dust 
medium. 

To understand the observed shape of the power spec- 
trum, we attempted to analyze the impact of variations 
in the beam size on the pattern of the observed power 
spectrum of the OMC1 image. We performed power 
spectrum calculations using three different aperture di- 
ameters: 18", 36", and 56". Here, the aperture diam- 
eter is an average for the observed region. The pixels 



in the data have been averaged using IMBIN task in 
MIRIAD (Sault et al. 1995). Each resolution of the 
newly derived map, 36" and 56" was achieved by av- 
eraging 9 and 14 pixel squares respectively. Each new 
pixel was replaced with the average value of the pixels 
within the resolution of the map. Despite differences 
in the details, the power spectra of the 1.1mm emis- 
sion at apertures of different sizes are similar in shape. 
However, as the size of aperture increases, the power 
decreases significantly (Fig. 2). The power spectrum is 
slightly less steep towards the large beam at relatively 
small scales, compared to the large scale structure of 
OMC1 at 1.1mm. The slope of the observed power 
spectrum is more affected after beam convolution than 
after averaging over the pixel values. For beam con- 
volutions, a gaussian beam was used to convolve the 
image at different resolutions using the CONVOL task 
in MIRIAD. The power spectra at lower and higher 
resolutions calculated from the beam convolution have 
similar appearances at spatial frequencies smaller than 
k w 30. But the slope changes significantly at higher 
spatial frequency as the aperture diameter increases. 
The average error in the noise was also measured in 
terms of the spatial power spectrum, and is presented 
in Fig. 2. The noise includes 1/f noise (Wilson et 
al. 2008) and confusion noise from the background. 
The shape of the power spectrum of the noise follows 
a power law and the present analysis yields a spectral 
index of 7 w —0.6. 

5. SUMMARY 

We present a 1.1mm emission map of the OMC1 re- 
gion observed with AzTEC, a new large- format array 
composed of 144 silicon-nitride micromesh bolometers 
that was in use at the JCMT. We performed a Fourier 
analysis of the image with a modified periodogram. 
The periodogram was measured by computing the dis- 
crete Fourier transformation. From the present anal- 
ysis of the OMC1 filaments at 1.1mm emission, the 
power spectrum steepens at relatively smaller scales. 
At larger scales, the power spectrum flattens and the 
large scale power law becomes shallower. The loga- 
rithmic slopes of the power-law spectra obtained seem 
to be lower than the logarithmic slope of Kolmogorov 
spectrum. We also performed power spectrum calcula- 
tions using three different aperture diameters. As the 
size of aperture increases, the power decreases signifi- 
cantly. We were able to fit at least three components of 
power law in the power spectrum of the 1.1mm emission 
map. The slope of the best fit at spatial frequency of 
30 £ k <; 100 is 7 « -2.66 ± 0.3. This result is similar 
to the spectral index of the power spectrum, 7 « —2.7 
that was found in numerical simulations. The effects 
of noise on the slope of the power spectrum were also 
included in the present analysis. 
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